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ABSTRACT 

We show that the two basic assumptions of the model recently proposed by Manrique and coworkers 
for the universal density profile of cold dark matter (CDM) halos, namely that these objects grow 
inside out in periods of smooth accretion and that their mass profile and its radial derivatives are 
all continuous functions, are both well understood in terms of the very nature of CDM. Those two 
assumptions allow one to derive the typical density profile of halos of a given mass from the accretion 
rate characteristic of the particular cosmology. This profile was shown by Manrique and coworkers 
to recover the results of numerical simulations. In the present paper, we investigate its behavior 
beyond the ranges covered by present-day A-body simulations. We find that the central asymptotic 
logarithmic slope depends crucially on the shape of the power spectrum of density perturbations: it is 
equal to a constant negative value for power-law spectra and has central cores for the standard CDM 
power spectrum. The predicted density profile in the CDM case is well fitted by the 3D Sersic profile 
over at least 10 decades in halo mass. The values of the Sersic parameters depend on the mass of 
the structure considered. A practical procedure is provided that allows one to infer the typical values 
of the best NFW or Sersic fitting law parameters for halos of any mass and redshift in any given 
standard CDM cosmology. 

Subject headings: cosmology: theory - dark matter - galaxies: halos 



1. INTRODUCTION 

The universal shape of the spherically averaged density 
profile of relaxed dark halos in high-resolution A-body 
simulations is considered one of the major predictions of 
standard cold dark matter (CDM) cosmologies. Down to 
1 % of the viri al radius R, it is wel l fitted by the so-called 
NFW profile (|Navarro et al.lll996t Il997t ). 



p(r) 



r(r s + r) 2 



(1) 



specified by only one mass-dependent parameter, the 
scale radius r s or equivalently the concentration c = 
r s /R. 

At radii smaller than 1 % of the virial radius, however, 
the behavior of the density profile is unknown. On the 
basis of recent numerical simulations, some authors ad- 
vocate a central asymptotic slope significantly steeper 
(Moore et al. 1998; Jing fc Suto 2000) or shallower 
(iTavlor fc Navarrdl200lt lRicottil[200l lHansen fe Stadell 
[2006D than that of the NFW law. Others suggest an 
ever decreasing abs o lute value of t h e logarithm i c slop e 
(|Power et all l200l E avarr 3 l200i iReed et alJ l2005h . 
which might tend to zero as a power of the radius as 
in the three-dimensional (3D ) Sersic (1968) or Einast o 
(Einasto & Haud 1989) laws (jMerritt et al.ll2005l[2006T) . 

This uncertainty is the consequence of the fact that 
the origin of such a universal profile is poorly under- 
stood. Two extreme points of view have been envis- 
aged. In one of these , it would be c aused b y repeated 
significant mergers (ISver fc Whitd 11998 : IRaig et al.l 



1998; ISubramanian et al.ll2000t iDekel et all [2003). while 



in the other it would be essentia lly the result of 
smoo t h accretion or se c ondar y infall (| Avila- Reese et all 
[1991 iNusser fc Shethl [19991 [DeTPopolo et all 1 2000 1 : 
Kulll Il999t iManriaue et al] 120031: IWilliams et all \20oi 
lAscasibar et alj|2004[ )~ 

The fact that the M — c relati on at z = is con- 
sistent (Salvador- Sole et al. 1998. IWechsler et all [2001 : 
IZhao et "all l2003bh with the idea that all halos emerge 
from major mergers with similar values of c, which then 
decreases according to the inside-out growth of halos dur- 
ing the subsequent accretion phase, seems to favor an im- 
portant role of mergers. However, the purely accretion- 
driven scenario is at least as attractive, as the inside-out 
growth during accretion leads to a typical density pro- 
file that appears roughly to have the NFW shape with 
the correct M — c relation in any epoch and cosmology 
analyzed (Manrique et al. 2003, hereafter M03). 

Certainly the effects of major mergers cannot be ne- 
glected in hierarchical cosmologies, so both major merg- 
ers and accretion should contribute in shaping relaxed 
halos. However, as noted by M03, if the density profile 
arising from a major merger were set by the boundary 
conditions imposed by current accretion, then the den- 
sity profile of halos would appear to be independent of 
their past aggregation history, so halos could be assumed 
to grow by pure accretion without any loss of generality. 
All the correlations shown by relaxed halos in numeri- 
cal simulations can be reco vered under this point of view 
([Salvador-Sole et al.l l2005h . Simultaneously, this would 
explain why halos with very different initial conditions 
and aggregation histories have similar density profiles 
([Romano-Diaz et aljl2006l ). 
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In the present paper, we show that the M03 model 
relies on two basic assumptions, namely, (1) that halos 
grow inside out in periods of smooth accretion, and (2) 
that the mass profile and all its derivatives are continuous 
functions. The former assumpti on is supported by the re- 
sults of numerical simulations (ISalyador-Sole et al.ll2005t 
ILu et al.ll2 006; Ro mano-Diaz et al.ll2006tl2007f) . while the 
second one is at least not in contradiction with them. In 
the present paper we show that both assumptions are 
in fact sound from a theoretical point of view, as they 
can be related to the very nature of CDM. This renders 
the predictions of the model beyond the range of current 
simulations worth examining in detail, as is done below. 

For simplicity, we consider spherical structures, which 
at best is an approximation to the triaxial structures 
observed in numerical simulations. Secondary infall is 
known to be influenced by deviations from spheric al sym- 
metry (|Bond fc Mv erd li996t IZaroubi et all 1 19961) . Yet, 



in an accompanying paper (fGonalez-Casado et al.ll2007l) . 
we show that this does not seem to affect the fundamen- 
tal role of the two assumptions given above. Another 
simplifying assumption used here is the neglect of sub- 
structure. It has recently been shown in high-resolution 
simulations that about 57 % of th e halo mass is col- 
lected in the previous major merger ( Faltcnb acher et al.l 
2005). However, substructures (and sub-substructures) 
contribute only abou t 5 % of the total mass in halos 
(jDiemand et alj [2007D whose density profile is well de- 
scribed by equation {]]). It seems therefore a good first 
approximation to ignore substructure. 

The paper is organized as follows. In £J21 we show how 
the M03 model emerges from the properties of standard 
CDM. The behavior of the predicted density profile at 
extremely small radii and for a wide range of halo masses 
is investigated in Our results are summarized in S2J 

2. CDM PROPERTIES AND HALO DENSITY PROFILE 
2.1. Inside-out growth during accretion 

We now argue that halos grow inside out during accre- 
tion, as is indeed found in numerical simulations, because 
of some properties of CDM; in particular, its character- 
istic power spectrum, leading to a slow halo accretion 
rate. As shown in £12.21 this has important consequences 
for the inner structure of these objects. 

Schematically, one can distinguish between minor and 
major mergers. In minor mergers, the relative mass in- 
crease produced, A = AM/M, is so small that the sys- 
tem is left essentially unaltered, whereas in major merg- 
ers A is large enough to cause rearrangements. 

The s maller the value of A, the most frequent the 
mergers (|Lacev fc ColelUQQl . For this reason, although 
individual minor mergers do not affect the structure of 
halos, their added contribution yields a smooth secular 
mass increase, the so-called accretion, with apparent ef- 
fects on the aggregation track M(t) of the halo. The 
accre tion-scaled rate, M/M(t), is given by (|Raig et al.l 
2001. hereafter RGS01) 



fta(M, t) 



dA A TZ m (M, t, A) , 



(2) 



where TZ m (M, t, A ) is the usual Lacey-Cole 
(|Lacev fc Cold |1993[ ) instantaneous merger rate (see 
eq. [Alj ^1 and A m is the maximum value of A for mergers 




Fig. 1. — M — c relation at z = predicted in the concordance 
model by the M03 model for A m = 0.26 (solid line) compared to 
the empirical rel ation traced by tho se points, with minimal error 
bars obtained by Zhao ct al. (2003a) from high-resolution simula- 
tions (triangles with error bars). 

contributing to accretion. In contrast, less frequent 
major mergers yield notable sudden mass increases or 
discontinuities in M(t). 

After undergoing a major merger (and virializing), ha- 
los evolve as relaxed systems until the next major merger. 
The fact that standard CDM is nondecaying and non- 
self-annihilating 1 guarantees that the mass collected dur- 
ing such periods is conserved. This does not yet imply 
that halos grow inside out during accretion, because their 
mass distribution might still vary due to energy gains or 
losses or to the action of accretion itself. Even if each 
individual minor merger would leave the halo unchanged 
their collective action might alter these systems. How- 
ever, standard CDM is also dissipationless and therefore 
halos cannot loose energy. Furthermore, under the as- 
sumption of spherical symmetry halos cannot suffer tidal 
torques from surrounding matter, and hence the sur- 
rounding matter is unable to alter the kinetic energy 
of the halo. Therefore, the only possibility for a time- 
varying inner mass distribution of accreting halos is that 
the accretion process causes it itself. 

This possibility, that the process of accretion itself 
could alter the internal mass distribution, would be re- 
alized only if the accretion time 1/1Z& were smaller than 
the dynamical time r cr . On the contrary, if is sub- 

stantially larger than r cr , the adiabatic invariance of the 
inner halo structure will be guaranteed, and the halo will 
evolve inside out. Thus, by requiring l/lZ a to be C times 
larger than r cr , we are led to the equation 



C 



-in 



GA vh (t)p(t) 



1/2 



K & (M,t), 



(3) 



which gives the upper mass M a for inside-out growth at 
t. In equation (J3]), p is the mean cosmic density and 

1 Both decay and annihilation rates are extremely small for re- 
alistic dark matter particle candidates. 
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Fig. 2. — Predicted density profiles (solid lines) compared to their fits to a NFW profile (dashed lines) for halo masses at z = 0, ranging 
from 10 16 Mq to 10 11 Mq (or, equivalently, from 10 3 to 10 — 2 times the current critical mass M*o for collapse), in the concordance model. 
Upper subpanels: Profiles. Lower subpanels: Residuals of the logarithmic fits. The halo mass and the corresponding best-fitting value of 
the NFW concentration parameter are quoted in each panel. 

A v j r (t) is the viria l izatio n density contrast given e.g. by 
iBrvan &: Normaiil (|1998l ). Here M a is indeed an upper 
limit because 1Z & (M, t) is an increasing function of t; see 
RGS01. In any CDM cosmology analyzed, equation 
appears to have no solution for values of C significantly 
larger than unity in the relevant redshift range. This 
means that accretion is always slow enough for halos to 
grow inside out, as required by the M03 model. 

As previously mentioned, the inside-out growth of ha- 
los in accreting periods is unambig uously confirmed by 
the r e sults of TV-body simulations (ISalvador-Sole et al.l 
120051 iLu et all l2006t iRomano-Diaz et al.1 120061 . 120071) . 
It is also consistent with the fact that dark matter 
structures preserve the memory of initial conditions, 
in the sense that the most initially overdense regions 
end u p being the central regions of the final struc- 
tures (Dicma nd et al.l [2005) , implying that the spatial 
positions of particles are not significantly perturbed by 
merging/accretion during the assembly of the structures. 
Likewise, the energy of the individual particles in the 
final structure (at z — 0) is very strongly correlated 
with their energies at much earlier times (z = 10; 
iDantas &; Ramoa l2006'l. This shows that particles even 
preserve the memory of the initial energies. 



2.2. Smoothness of the mass profile 

Contrarily to an ordinary fluid, CDM is collisionless 
and free- streaming and, hence, cannot support disconti- 
nuities (shock fronts) in the spatial distribution of any of 
its macroscopic properties. As a consequence, all radial 
profiles in relaxed halos are necessarily smooth. This 
holds in particular for the mass profile, M(r), and its 
radial derivatives, 2 which has the following consequence. 

The inside-out growth of a halo during accretion f ^2.ip 
implies that the mass profile M(r) built at that interval 
is the simple conversion, through the definition of the 
instantaneous virial radius 

1/3 



R(t) = 



3M(t) 



(4) 



AirA vil (t)p(t) 

of the associated mass aggregation track M(t). 

The smoothness condition implies that the old M(r) 
profile must match perfectly with the new part of the 
profile built during that time. Since minor mergers only 
cause tiny discontinuities, the new piece of the M(t) track 
that they produce is well approximated by a smooth 

2 We are naturally referring to the theoretical mass profile, equal 
to the average profile over random realizations or over time real- 
izations as long as the halo evolves inside-out. 
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function, and, as the functions p(t) and A v i r (t) in equa- 
tion Q are also smooth functions, the corresponding 
piece of the M(r) profile automatically fulfills the right 
smoothness condition. Thus, the system can grow, dur- 
ing accretion, without the need to essentially rearrange 
its structure. 

Only when a halo undergoes a major merger and its 
M(t) track suffers a marked discontinuity will the mass 
profile prior to the major merger no longer match the 
piece that begins to develop after it. Since the M(r) 
profile cannot have any discontinuities, the halo is then 
forced to rearrange its mass distribution (through violent 
relaxation) to fulfill the required smooth condition. 

In other words, the fundamental assumption of the 
M03 model that the mass distribution of relaxed halos 
is determined by their current accretion rate (through 
dramatic rearrangements of the structure on the occa- 
sion of major mergers and very tiny and negligible ones 
during accretion periods), is simply the natural conse- 
quence of the slowly accreting, nondecaying, nonself- 
annihilating, dissipationless, and collisionless nature of 
standard CDM. 

2.3. The M03 model 

The mass profile of a specific halo with mass M\ at 
time ti accreting at a given rate during any arbitrar- 
ily small time interval At around t\ is therefore simply 
the smooth extension inward of the small piece of pro- 
file being built during that interval. 3 Unfortunately, the 
smooth extension of a small piece of a function is hard 
to find in practice, so the mass profiles of real individual 
halos can hardly be obtained in this way. 

There is one case, however, in which such a smooth 
extension can readily be achieved: that of halos with Mi 
at t{ accreting at the typical cosmological ratelZ a (M(t),t) 
during any arbitrarily small interval of time around ti . In 
this case, the (unique) smooth extension we arc looking 
for necessarily coincides with the smooth function M(t), 
the solution of the differential equation 

7l a (M(t),t) 



M{t) 



(5) 



for the boundary condition M(t\) = Mi, properly con- 
verted from t to r by means of equation (j4]). Once the 
typical M{r) profile is known, by differentiating it and 
taking into account equations J!]) and ((5|, one is led to 
the typical density profile for halos with Mi at ti proposed 
by M03, 



p{r) = 



5{t) = 




(6) 



(7) 



K a {M{t),t) dt 

From equations © and ([7]) we see that the shape of 
this profile is ultimately set by the CDM power spec- 
trum of density perturbations in the cosmology consid- 
ered through the merger rate lZ m (M, t, A) used to calcu- 
late the accretion rate 7?. a (M, t) (see eqs. |A1| and [2]). 



3 The derivatives at any order of a smooth function, which is 
known in some finite domain, are completely determined at any 
point of that domain. Thus, by taking the Taylor series expansion 
of the function in that point, one can extend it in a unique way 
outside the domain. 
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Fig. 3. — Relations between the NFW r B and M s parameters 
obtained in the concordance model [sghown with solid lines for 
masses up to 10M* (z) and dashed lines beyond it] from the fits of 
the theoretical density profiles of halos with varying mass at the 
redshifts 0.0, 1.03, 2.08, and 3.94. The dot-dashed line gives the 
analytical fit (see text) to the mean of the solid lines. 

This dependence is, however, so convoluted that the den- 
sity profile (eq. [6]) must be inferred numerically. Only 
its central asymptotic behavior can be derived analyti- 
cally, as will be shown in the next section. 

3. SOME CONSEQUENCES OF THE MODEL 

From equation © we see that the exact shape of the 
density profile in equation ([6]) depends on A m . This pa- 
rameter marks the effective transition between minor and 
major mergers, and it can be determined from the em- 
pirical M — c relation at some given redshift. 

For each given A m value, the density profiles, down 
to -R/100, predicted at z = in the concordance model 
characterized by (f2 m , f^A, h, cr s ) = (0.3,0.7,0.7,0.9) for 
halos with different masses have been fitted to the NFW 
profile to find the best-fit values of c. Then we searched 
for the value of A m that minimizes the departure of the 
theoretical M — c relations from the empirical one draw n 
from high-resolution simulations by I Zhao et al.l (|2003a). 
As shown in Figure [TJ A m = 0.26 gives an excellent fit 
over almost 4 decades in mass (8 x 10 10 h~ 1 M & < M < 
4 x lO^h' 1 M Q ). 

In Figure [21 we plot, down to a radius equal to the cur- 
rent resolution radius of most numerical simulations, the 
density profiles predicted in the concordance model for 
halo masses at z — ranging from 10 11 to 10 16 M Q . They 
are all well fitted to a NFW profile, although there is a 
tendency for the theoretical profiles for M > 10 14 M to 
deviate from that shape and approach a power law with 
logarithmic slope intermediate between the NFW asymp- 
totic values of —1 and —3. This tendency also makes c 
increase very rapidly at large masses where the fit by the 
NFW law is no longer acceptable. This causes the M ~c 
relation to deviate from its regular trend at smaller val- 
ues of M (see Fig. 2). Both effects, already reported in 
M03, were later observed in simulated halos ; Zhao et al.l 



DENSITY PROFILE AND CENTRAL BEHAVIOR OF HALOS 



1^ 

\ 

3 

O 



< 



la 

O 



< 



12 
9 
6 
3 

10 

10 

12 
9 
6 
3 

10 

10 



TTTTTTTTTTTTTTTTTTT 

c=1.57xl0 24 - 
n=15.7 : 



-M=10 16 M ( 

I I I I I I I I I I I I I I I I | : 



: l I I I I I I I I I I I I I I I : 



c =1.79xl0 7 " 
n = 5.69 




I I I I I I I I I I I I I I I | : 



TTTTTTTTTTTTTTTTTTT 

c =1.42xl0 15 - 



n=10.9 




-M=10 14 M 

1 I I I I I I I I I I I I I I I I I I 




: l I I | I I I I I | I I I I I | I I I : 

c =1.52xl0 6 ~ 
n=4.87 : 



^M=10 B M, 

- I I I I I I I I I I I H-H+H-H 



i i i I I Im iti ill I Im i i I i i i i i I i i i i i I i i i 



■i i i | i i i i i | i i i i i | i i i- 

c =3.78xl0 9 ~ 
n=7.34 




I I I I I I I I I I I I I I I | : 



: l I I | I I I I I | I I I I I | I I I : 

c =5.30xl0 5 ~ 
n=4.52 : 



^M=10 6 M, 

- I I I I I I I I I I I hH+H+H 



0.3 0.6 0.9 



0.3 0.6 0.9 
(r/R) 1 /" 



0.3 0.6 0.9 



Fig. 4. — Same as Figure [2] but for the Sersic profile. All the profiles are now drawn down to a 1 pc radius and cover a wider range of 
halo masses: 10 16 - M Q (or, equivalently, from 10 3 M* to 10~ 7 M* )- 



20033 iTasitsiomi et al.ll2004[) . This is a clear indication 
that the NFW profile does not provide an optimal fit for 
very massive structures. Of course, above 10 14 M Q halos 
are hardly in virial equilibrium, so such a deviation has 

essentially no pra ctical effects. 

As explained in [Salvador-Sole et alj (|2005f ). another in- 
teresting consequence of the M03 model is that the M(t) 
tracks traced by accreting halos (hence, growing inside 
out) coincide with curves of constant r s - and M s - val- 
ues, with M s defined as the mass interior to r a . Thus, 
the intersection of those accretion tracks at any arbitrary 
redshift sets the relation M s (r s ) between such a couple 
of parameters, implying that the M s (r s ) relation satis- 
fied by halos is time-invariant. In Figure[31 we show how 
the different M s (r s ) curves obtained by fitting the den- 
sity profiles predicted at different redshifts to a NFW law 
overlap. There is only some deviation at large masses, 
where the density profiles are not correctly described by 
the NFW profile. Such a time- invariant M s (r s ) relation 
is well fitted, for r s in the range 10~ 4 Mpc < r s < 1 Mpc, 

by 



Substituting equation ([5| into the relation 

M s _ ln2-0.5 
IT ~ g(c) 

which holds for NFW profiles, one is led to 



47rA vir (£) p(t)~ 


0.82 


M 


0.18 / s, 

5(c) 


2.37 x 10 6 p _ 




1O 13 M _ 


c 2.45 



Ms 



10 13 M, 



141 







Mpc 



2.45 



(8) 



(9) 



(10) 



where g(c) stands for ln(l + c) — c/(l + c) and po is 
the current mean cosmic density. Equation (|10j) is an 
implicit equation for the concentration of halos with any 
given mass and redshift. 

What about the central behavior of the predicted den- 
sity profile? According to equation (H|) , small radii corre- 
spond to small cosmic times. In this asymptotic regime, 
all Friedman cosmologies approach the Einstein— de Sit- 
ter model in which A V i r (t) is constant and pit) is (in the 
matter-dominated era when halos form) proportional to 
t~ 2 . If the power spectrum of density perturbations were 
of the power-law form P{k) oc fc J , with the index j sat- 
isfying 1 > j > —3 to guarantee hierarchical clustering, 
then the universe would be self-similar. The mass ac- 
cretion in equation ([5]) would take the asymptotic form: 

M(t) cx tlta ( S ee the Appendix). The fact that both 
A V n (t)p(t) and M(t) would then be power laws has two 
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consequences. First, the time dependence of their re- 
spective logarithmic derivatives on the right-hand side 
of equation ([7]) cancel, which implies that p(t) is pro- 
portional to A v i T (t)p(t), and hence to p(t). Second, the 
virial radius given by equation (j4|) is also a power law, 
R(t) cx [M(t)/p(t)]V3 cx t 2 a+ 4 )/[3W+3)l. From this we 
get t(r), and thereby one finds 



3(3 + 3) 

p(r) cx r 5+3 



(11) 

This central behavior, fully in agreement with the nu- 
merical profiles obtained from power-law power spec- 
tra, is particularly robust, as it does not depend on 
A m . Note that it coincides with the solution derived 
in self-similar cosmologies bv lHoffman fc Shahaml (|1985f ) 
assuming spherical collapse. It is worth noting, however, 
that in that derivation, such an asymptotic behavior is 
restricted to j > — 1 so as to warrant the req uired adi- 
abatic invariance (jFillmore fe Goldreich| [T984). while, in 
the present derivation, there is no such a restriction as 
the central density profile is not assumed to be built by 
spherical infall, but results instead from smooth adapta- 
tion to the boundary condition imposed by current ac- 
cretion. 

The CDM power spectrum is, of course, not a power 
law. However, in the limit of small masses involved in 
that asymptotic regime, it tends to a power law of index 
j = —3. Thus, according to equation ([LI]) , we expect a 
vanishing central logarithmic slope of p{r) for the stan- 
dard CDM case. This is confirmed by the numerical pro- 
files obtained in this case: as one goes deeper and deeper 
into the halo center, they become increasingly shallower. 

What is still more remarkable is that down to a radius 
as small as 1 pc, the density profiles appear to be well 
fitted by the 3D Sersic or Einasto law 



p(r) = p exp 



l/ri 



(12) 



over at least 10 decades in halo mass (see Fig. [4]), from 
10 6 M Q to 10 16 M Q . We remind the reader that for 
power-law spectra, equation © leads to density profiles 
with central cusps, so the Sersic shape is not a general 
consequence of the M03 model, but it is specific to the 
standard CDM power spectrum. In fact, from the rea- 
soning above we see that what causes the zero central log- 
arithmic slope in the CDM case is the fact that the log- 
arithmic accretion rate dlnM(£)/dlni = t1Z a (M(t),t) 
diverges in the limit of small values of t. This is in con- 
trast to the general power-law case, where the accretion 
rate remains finite. 

Similar to the characteristic density p c in the NFW 
profile (eq. [Tj), the central density po entering the 3D 
Sersic law (eq. |12j ) can be written in terms of the mass 
M and the values of the two (instead of one) remaining 
parameters, n and either r n or c n = R/r n : 



M 



PQ = 



47mr n L(3n) 
where T is the usual gamma function and 



P{a,x) = 



1 



df e"* t"-- 1 



(13) 



(14) 



r(a) jo 

is the incomplete or regularized one. At z = 0, the 
two free parameters n and c n depend on M according 




10 12 14 
log(M/M Q ) 



16 



Fig. 5. — Mass dependence at z = of the Sersic law parameters 
n and c n fitting the predicted density profiles of halos in the con- 
cordance model (solid lines) and their fits by the simple expressions 
given in the text (dashed lines). 

to the relations plotted in Figure which are well ap- 
proximated by 



n(M) = 4.32 + 7.5 x 10" 



ln[c n (M)] = 13.3 + 7.5 x 10" 



In 



M 



In 



M 

~Mr, 



4.6 



n 5.6 



leading to the following combined relation 



c n (M) = 5.84 x 10 5 



M 



(15) 



(16) 



(17) 



These expressions can be used to infer the typical values 
of the Sersic parameters for present-day halos with any 
mass. It is worth mentioning that the predicted values of 
n are of the same order of magnitude as the ones obtained 
bv lMerritt et al.l (|2005l ) from simulated halos with masses 
ranging from dwarf galaxies to galaxy clusters (the values 
of c n were not presented in that work). 

To obtain more general values of these parameters for 
halos of any mass and redshift, we can proceed as in the 
NFW case above. For reasons identical to the ones lead- 
ing to the time- invariant relation M s (r s ), the relations 
Po(r n ) and n(r n ) must be time- invariant. In Figure O 
we see how the corresponding curves obtained from the 
fit to the Sersic profile of the same density profiles as 
used in Figure [3] overlap, indeed, even better than the 
M s (r s ) curves do, since there are no large deviations at 
large masses. These invariant relations are well fitted, 
for x = log(r n /Mpc) in the range —25 < x < —9, by 



Po(r n ) = Poexp(A) , 
A = -2.48- 4.73x 



(18) 



-0.270x 2 



6.24 x 10~ 3 x 3 . 
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-24 -21 -18 -15 
log(r n /Mpc) 

Fig. 6. — Same as Figure[3]but for the po( r n) and n(r n ) relations 
among the best-fit Sersic parameters. The respective analytical fits 
(see text) to the mean curve are shown in dot-dashed lines. Here 
po is the current mean cosmic density. 



n(r n ) = -4.43 - 1.43a; 



(19) 



-Q.0313ar - 2.98 x 10 



Replacing these expressions into equation (|13|) . we can 
solve for r n and then use the equation IT9|) to find the 
value of n. 

This provides a very concrete prediction that can be 
tested with numerical simulations. One can take the very 
strong correlation shown in Figure [6l which allows one to 
fit the density profile of any dark matter structure with 
only two free parameters: e.g. n and pq. With this value 
of n (purely from the shape of the density profile), one 
now gets a value for the mass (from Fig. [5]) . This mass 
can trivially be compared to the true virial mass (which 
is naturally known in the simulation), and hence one can 
confirm or reject the prediction of the accretion-driven 
model. 



shown to be in overall agreement with the results of nu- 
merical simulations. In this model, the density profile of 
relaxed halos permanently adapts to the profile currently 
building up through accretion and does not depend on 
their past aggregation history. As a consequence, the 
typical density profile of halos of a given mass at a given 
epoch is set by their time-evolving cosmology-dependent 
typical accretion rates. 

Although halos have been assumed to be spherically 
symmetric throughout the present paper, this is not cru- 
cial for the M03 model. As will be shown in a follow- 
ing paper (Gonzalez-Casado et al. 2007), the results pre- 
sented here also hold for more realistic triaxial rotat- 
ing halos. Furthermore, an approach similar to the one 
followed here allows one to explain not only their mass 
distribution, but also others of their structural and kine- 
matic properties, such as the radial dependence of angu- 
lar momentum. 

According to the M03 model, the central asymptotic 
behavior of the halo density profile depends, through the 
typical accretion rate, on the power spectrum of den- 
sity perturbations. The prediction made in the case of 
power-law spectra should be possible to check by means 
of numerical simulations, provided one concentrates on 
massive halos, as these reach the asymptotic regime at 
larger radii. In the case of the standard CDM power spec- 
trum, the model predicts a vanishing central logarithmic 
slope. The way this asymptotic behavior is reached is 
surprisingly simple: down to a radius as small as 1 pc, the 
density profile is well fitted by the 3D Sersic or Einasto 
profile over at least 10 decades in halo mass. 

Another consequence of the M03 model with useful 
practical applications is the existence of time-invariant 
relations among the NFW or 3D Sersic law parameters 
(M s and r s in the former case and po, r n , and n in the 
latter) fitting the halo density profiles. A code is publicly 
available 4 that computes such invariant relations for any 
desired standard CDM cosmology. 

Some of these consequences can be readily tested by 
numerical simulations or by (X-ray or strong- lensing) ob- 
servations, which should allow one to confirm or reject 
the prediction for the central behavior of the density pro- 
file of halos that is made by the M03 model. 



4. SUMMARY 

We have shown how the basic properties of standard 
CDM can justify the M03 model, which was previously 
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AYA2006-15492-C03-03. We thank Donghai Zhao and 
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1 http : //www . am . ub . es/cosmo/gravitat ion . html 
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APPENDIX 

THE ACCRETION TRACK ASYMPTOTIC BEHAVIOR 



The instantaneous rate of mergers undergone by halos with mass M at time t yielding a relative mass increase of 
A, is (lLacev fc Coldfl993h 



K m (M,t,A) 



y/2/nM 
a 2 {M') 



d5 c da(M') 



dt dM> 



1 



a 2 (M') 



a 2 (M) 



-3/2 



exp 



s 2 c (t) 

2<j 2 {M') 



1 - 



a 2 (M>) 



a 2 (M) 



(Al) 



where M' = M(l + A) is the final mass, S c (t) is the linearly extrapolated, cosmology-dependent, critical density 
contrast of density fluctuations at t for collapse at the present time to, and a(M) = cr(M, to) is the rms density 
contrast of fluctuations on scale M at to, related to the zeroth-order spectral moment. 
In a self-similar universe with spectral index j, one has 

.3+3 



<\ A. I ) = <\ , , ( ^ 



-2/3 



a(M) = a(M+ ) 



M 
M^o~ 



(A2) 



where S c q and M+q are the current values of S c (t) and Af*(i), respectively, M±(t) being the critical mass for collapse 
at t, the solution of the implicit equation a(M+,t) = S c q. Therefore, equation (|A1[) takes the form 



K m (M,t, A) 



2 j + 3 v{M,t){l + Af^r- 

7T 9* f,. .,,+3 13/2 



v(M, t) = 



(1 + A) 

tS -2/3 

ta 



■ exp 



- 1 
M 



v 2 (M, t) 



(1 + A) 



M 



M*(t) 



(A3) 



(A4) 



Substituting equation (|A3|) into equation ^ and changing A by x = ^ 2 [(1 + A)^ +3 '/ 3 — 1], we obtain the following 
expression for the accretion rate 



K a (M,t) = 



ZAv 2 
(7+3)* Jo 



j(M,i,A m ) 



d.r 



{l + v- 2 x)^ 



, i+ expj-z/2) 



where, for simplicity we have dropped the explicit dependence of v and used the notation A = \/2/ir (J + 3)/9 and 
x m (M,t,A m )=v 2 [(l + A m )V+ 3 V 3 -l]. 

In the differential equation (eq. [5]), the variable M in lZ a is replaced by the mass accretion track M(t). As halos 
grow through both accretion and major mergers, the accretion tracks M{t) increase with increasing time less rapidly 
than do M*(t), tracing the typi cal m ass evolution of halos in any self-similar universe. Thus, v{t) = v{M(t),i) is a 
decreasing function of t (see eq. |A4| ) and, in the small-* asymptotic regime, v{t)~ 2 tends to zero. Taking the Taylor 
series expansion of [1 + i/ _2 (t)x] 3 /^ +3 - ) inside the integral on the right-hand side of equation (|A5|) . at leading order in 



L (*), we obtain 



K a (M(t),t) = 



9AV 



(J + 3) 2 



'"f" , exp(-x/2) 
dx = 



9V2^Ar 
U+3? 



erf 



Ktwa+Am)^ 



(A5) 
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where eri(x) is the error function. Given the value of constant A, and given that, for v{t) tending to infinity, the error 
function approaches unity, we are led to 

n a {M{t),t) = ^-r 1 . (A6) 

J + 3 

Finally, integrating equation J5J for such an accretion rate, we are led to the following asymptotic behavior for accretion 
tracks: 

M(t) oc t~ . (A7) 

4 

Note how it compares with the time dependence of the critical mass in a self-similar universe: M+(t) oc ti+ 3 . 



